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Abstract 

We examine the one-dimensional spinless Falicov-Kimball model extended by a hybridization 
potential between the localized and itinerant electron states. Below half-filling we find a crossover 
from a mixed-valence metal to an integer-valence phase separated state with increasing on-site 
Coulomb repulsion. This crossover regime is characterized by local competition between the strong- 
and weak-coupling behaviour, manifested by the formation of an electronic polaron liquid. We 
identify this intermediate-coupling regime as a charge-analogy of the Griffiths phase; a phase 
diagram is presented and discussed in detail. 

PACS numbers: 71.28.+d, 71.27.+a 
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The Falicov-Kimball model (FKM) has found widespread use as a model of the charge 
physics in mixed-valence (MV) materials. Although the FKM's description of a band of 
itinerant conduction (c-)electrons interacting via a repulsive contact potential G with an 
underlying lattice of localized (/-)electrons retains its relevance,— the great evolution of 
our understanding of MV systems has required modification of the model's original form.— 
The most popular such modification is the assumption of some overlap between the c- and 
/-wavefunctions: this yields the extended or quantum Falicov-Kimball model (QFKM). 

The ground state of the QFKM has for many years been a matter of controversy, with 
two contradictory pictures common in the literature. The homogeneous solution 3 assumes 
a simple renormalization of the G = band structure by the Coulomb interaction; at mean- 
field level, this has been used to deduce the existence of a novel electronic ferroelectric MV 
state.— On the other hand, the relationship of the single-impurity QFKM to the Kondo 
impurity 5 suggests distinct strong- and weak-coupling regimes in the lattice models 

The latter scenario leads to interesting properties of the intermediate-coupling regime, 
specifically the formation of electronic polarons.^ This electronic polaron model has re- 
cently been used to explain the transport properties of certain Heusler and Ce alloys?^ more 
intriguingly, such polaronic structures play a pivotal role in the theory of the exotic "valence- 
fluctuating" superconductivity in CeCu2Ge2 and CeCu2Si2.— It is therefore of considerable 
interest to study the formation of electronic polarons in the QFKM. 

In this letter we address this question by a non-perturbative study of the one-dimensional 
(ID) QFKM below half- filling. We derive an effective Hamiltonian for the occupation of the 
/-orbitals which gives a complete description of the physics. This allows us to exactly study 
the intermediate-coupling regime, where we find local competition between the strong- and 
weak-coupling phases. This produces a charge analogy of the Griffiths phase, which may be 
visualized as a gas of electronic polarons moving in a MV background. A schematic phase 
diagram is presented and discussed in detail. 

The ID QFKM for spinless fermions has the Hamiltonian 

n QFKM = E 4c,- + £f E 4 + G E n H + { c \h + H - c -} W 

(hj) J 3 3 

Some overlap between the c- and /-electron wavefunctions is assumed hence the on-site 
hybridization term V -C t. The concentration of electrons is fixed at n = (1/N) J2j {j 1 ^ + n j} 
where N is the number of sites. In general we have n < 1 due to the presence of doped 
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impurities. We assume a MV state at G = 0, requiring the Fermi level to be pinned at 
— 2t < 6f < — 2tcos(n7r) in the noninteracting system. 2 

Bosonization is an asymptotically exact method, rigorously describing the long- 
wavelength collective behaviour of the itinerant electrons;— the short-range behaviour re- 
mains fermionic, however, and a cut-off a on the wavelength of the bosonic density fluctu- 
ations is therefore imposed.— In simple cases such as the Hubbard model, a can be taken 
as vanishing. This is not true for systems such as the QFKM where the interactions with 
the localized orbitals determine the lattice constant as the minimum physical length scale. 
We are thus required to assume finite a > a: in the QFKM a parameterizes the electron 
derealization length. 

A finite wavelength cut-off may be easily included into the standard bosonization formal- 
ism.— In particular, for the linearized fermionic fields about the Fermi points the well-known 
Mandelstram representation still holds 

l~Aa - 

c uj = d—F u exp {-iv [<j)(xj) - u9(xj)\) (2) 

where and 6 are the dual Bose fields, and the subscript v = L(—), R(+) as subscript 
(otherwise) for the left- and right-moving fermions respectively. The Bose fields are defined 
in terms of particle-hole excitations about the two Fermi points. For a system of length 
L>awe have 

4>{*i) = -iEElf^MA)^ (3) 

o(*j) = -*££^p^)a^K^ (4) 

v k^O KLj 

where p u {k) = J2k' c uk'-k c v,k' are the particle-hole density operators. The wavelength cut-off 
is enforced in Eq. (jSJ) and Eq. (J3J) by the function A a (k): this function satisfies the conditions 
A a (k) pa 1 for \k\ < a" 1 and A Q ,(/c) w otherwise. The numerical constant A in Eq. (J2J is 
determined by the functional form of A a (k). F v is the so-called Klein factor . 11 ' 12 

Eq. (J2J) correctly reproduces the long-wavelength (^> a) fermionic anticommutators and 
correlation functions. Since the Bose fields cannot resolve separations less than a, however, 
the Mandelstram identity breaks down at short distances ~ a. This is reflected in the Bose 
field commutators, which are "smeared" by the cut-off function. For example, assuming 
exponential cut-off A a (k) = exp(— a|fe|/2) we have [4>(xj),9(0)]_ = —iaxctan(xj/a) and 
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[d x (p(xj),8(0)]_ = —ia[a 2 + (re,) 2 ] -1 . In the limit a — > we recover the Luttinger model 
forms of a sign and Dirac delta function respectively.— 

Substituting Eq. (J2J) into Eq. (JTJ) and applying standard field-theory techniques, it is a 
simple matter to obtain the bosonized form of the QFKM Hamiltonian. For insight into the 
physics, however, we must rotate the Hilbert space to couple the c- and /-electron orbitals. 
To accomplish this, we apply a lattice generalization of the shift transformation used by 
Schotte and Schotte in the X-ray edge problem, U = exp {i^^J2j( n j ~ ^^(^j)}-" The 
Fermi velocity is defined Vp = 2tasm(kpa) where kp = n^n/a, with Hq the noninteracting 
c-electron concentration. 

In the transformed Hamiltonian we combine the /-electron operators and the Klein factors 



into pseudospins using a generalized Jordan- Wigner transformation r 2 



± r+ 

2> '] 



f}F v e 



-il/TTXj /2q 



a exp (— i^j- sgn(xj' — Xj) (n^ — |U. Note the relationship of r 2 to the 
occupation of the /-orbital at site j. We then write the transformed Hamiltonian in its final 
form 

Vpd 



frn FKM u = ^-E{(^(^)) 2 + (^(^)) 2 } + g(^-|)E 
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+ H.<i) 



where K. a {xj) = S a (xj) + C a (xj), S a {xj) = ££°=i[{^ [0(x j+n ), (f>{xj)}_} - f] (r/ +n - r/_ n+1 ) 
and £ Q (^-) = - l) E~=i{^ (r/ +n - r/_ n ). We have kept all terms 

produced by the canonical transform in Eq. (jSJ. 

The three terms in Eq. (jSj) arising from the Coulomb interaction are of special note. The 
first term represents a renormalization of the /-level energy and drives a valence transition 
with increasing G. For < Uq < 0.5 (0.5 < Uq < 1) the negative (positive) sign of this 
energy-shift implies the emptying of the c-electron (/-electron) band so that all electrons 
have an unambiguously /-electron (c-electron) character. The case < < 0.5 is partic- 
ularly interesting as here the Ising interaction term orders the available /-electrons into a 
single contiguous block, the well-known segregated (SEG) phase of the V = limit. 15 The 
origin of this interaction is the finite spread a > a of the c-electron wavef unctions: to mini- 



mize the interorbital Coulomb repulsion this favours empty underlying /-orbitals. The last 
term originates from the backscattering of the c-electrons off the localized orbitals. In the 
absence of the hybridization, this will order the available /-electrons into crystalline phases 
via a Peierls-like mechanism.— The addition of the hybridization, however, replaces these 
phases by a MV state; in the following analysis we therefore neglect the backscattering. 
Including this effect does not qualitatively alter our conclusions. 17 

The phases of the QFKM are most usefully classified in terms of the occupation of the lo- 
calized orbitals. We therefore derive an effective Hamiltonian only for the /-occupation (i.e. 
the r-pseudospins) . This is straightforwardly achieved by replacing the Bose fields in Eq. (JHJ) 
with suitably chosen expectation values. At weak-coupling, Schlottmann's renormalization 
group study reveals that the system flows to the G = fixed point;— at strong-coupling, 
where the system is in the SEG phase, exact diagonalization calculations by Farkasovsky 
in the V = limit indicate that the c-electrons are at their noninteracting fixed point.— 
On the basis of these studies it is natural to take (<j)(xj)) = (9(xj)) = across the phase 
diagram. We hence obtain the effective Hamiltonian 

n ca = -JJ2 t*t; +1 + G (n* -§)$>; + T J cos (M*i) + W + f 1 xj) (6) 

3 j 3 ■ 

We have approximated the segregating interaction by its nearest-neighbour form where J = 
f lE fc A a (fc)e itfl . 

7i cB is instantly recognizable as a transverse field Ising model. For large G, we recover 
the SEG phase as the Ising term dominates the Hamiltonian. In the limit of vanishing 
Coulomb repulsion, however, the site-dependent transverse field h® = 4V cos(JC a (xj) + 
[kjF + dominates Eq. (jHJ). The pseudospins are therefore ordered along the x-axis 

with (tJ) ~ |sgn(/ij). Translated back into the /-occupation 'language' this implies that 
(fjCj + H.c.) 7^ 0, with the occupation of each /-site being given by (nj) w nl where Uq is 
the /-electron concentration for G = 0. This clearly corresponds to a MV state. 

The MV and SEG phases are characterized by quite contradictory behaviours of the /- 
electrons: in the former the localized orbitals are in a superposition of their different valence 
states, while in the latter there is a phase separation between regions with different integer 
ionic valence. At coupling strengths intermediate between the MV and SEG phases, there- 
fore, we expect to observe some competition between these two distinct charge physics. For 
a vanishing longitudinal field, Eq. (0) displays an order-disorder transition; the longitudinal 



field however lifts the system from criticality and we instead find a crossover between the 
ordered and disordered phases.— To understand the details of the crossover regime, we must 
examine the site-dependence of h x y, specifically, we study the operator fC a (xj), which we 
replace by its expectation value. 

The string operator JC a (xj) consists of two contributions, one measuring the short-range 
[«S a (xj)] and the other the long-range [C a (xj)} ordering of the pseudospins. In the SEG phase, 
the short-range correlations of the pseudospins are very strong and so S a (xj) vanishes; the 
importance of disorder in the crossover regime, however, means that here S a (xj) has a ran- 
dom variation. In contrast, the long-range ordering of the pseudospins in the SEG phase per- 
sists deep into the crossover regime. C a (xj) is only non-zero for < Uq < 0.5 where it takes a 
maximum at the edge of the /-electron block, increasing linearly as this point is approached 
from either side: for the boundary of the /-block at xq we find C a (xj) ~ (J^; ~ l) The 
linear variation of C a (xj) implies a quasiperidoic transverse field in the effective Hamilto- 
nian Eq. The periodicity of this field is in general incommensurate with the lattice: 
numerical 19 and analytical 20 studies have revealed the physics to be qualitatively identical 
to a random variation of transverse field. Fisher has used a real-space renormalization group 
method to exmaine in detail this random transverse field Ising model.— To utilise his de- 
tailed results we replace the transverse field in Eq. (JHJ) by a field which varies randomly 



with j, the values being drawn from a cosine distribution p{h)dh = (CV) Vl - (h/C) 



where C = AV^jAa/a. 

Following Fisher we find it convenient to define the dimensionless parameter 5 = 
\og(G c /G) where G 2 C = kt sin(riQ7r)V. Assuming exponential cut-off, the proportionality 
constant takes the form k = f + ^ • Although the dependence upon a obviously 

restricts the predictive powers of our analysis, it does not however prevent us from drawing a 
schematic phase diagram: for this purpose it is acceptable to choose a = constant. Without 
loss of generality we therefore set k equal to unity, being equivalent to renormalizing the 
Coulomb interaction G — > G/yk. In the QFKM's phase diagram presented in Fig. $IJ we 
plot the ratio G 2 /tV as a function of n^. We note that the phase diagram for 0.5 < Uq < 1 
can be obtained by reflecting Fig. (0) along the line = 0.5. The only significant dif- 
ference to Fig. (JTJ) is that the SEG phase is replaced by the "empty" state (tj/ = 0). We 
obtain the V — > weak-coupling crystalline phases in the phase diagram on the inclusion of 
the backscattering corrections in Eq. (j3J). Physically, since a finite hybridization is always 



6 



present, a MV state will be realized instead of these phases, as we have observed above. 

Starting in the segregated phase (5 <C 0) we lower G. Below the value G = ■Jtt/2G c , 
the Ising coupling in Eq. (jSJ) is no longer greater than the transverse field everywhere on 
the lattice. The rare regions where the transverse field term is strongest breaks the single 
/-electron block up into randomly distributed large clusters separated by regions of the 
MV phase: this local competition between the MV and SEG phases may be regarded as 
a charge analogy of the Griffiths phase [Gl in Fig. CQ)]. 21 Although such a state has been 
inferred by numerical work, this is the first analytic indication of its existence.— The system 
still, however, retains recognizable features of the SEG phase: the mean / — / correlations 
(njnj +x ) decay exponentially towards (n^) 2 > (nl) 2 with a correlation length £ ~ e~ 2S2 , 
the non-universal power-law dependence reflecting the continued presence of the infinite 
/-cluster throughout Gl. 

As G is lowered the Griffiths phase Gl gives way to the quantum critical regime (QCR). 
The QCR is dominated by the physics of the critical point: at every length scale we find the 
coexistence of equal regions of mixed- and integer- valence. Due to the /-level renormalization 
we do not find algebraic decay of correlation functions, but rather exponential attenuation 
with a correlation length of £ ~ 1 + 0(5 2 ). For = 0.5 the /-level shifting vanishes and so 
the QCR shrinks down to a single quantum critical point (QCP) characterized by dynamical 
critical exponent z = oo at G = G c (tIq = 0.5) = y/tV . Note that the precise position of the 
QCP is dependent upon the choice of <x— 

Further decreasing G below the QCR, we reach a second Griffiths phase (G2) with oppo- 
site character to Gl: here the hybridization term dominates the segregating interaction over 
most of the lattice. The SEG phase almost completely disappears, with only rare clusters of 
valence-ordered localized orbitals embedded in an MV background remaining of this state. 
The mean correlation functions decay exponentially, although less rapidly than in the Gl, 
with £ ~ (l+5~ 2 ) -1 . Since the minimum magnitude that the transverse field in 7i cS assumes 
is min \h®\ =0, the weakly disordered phase should be present for any G > 0; for 5 ^> 1 the 
integer- valence clusters are however extremely rare and the last two terms of 7i eS dictate the 
physics. As such, we can justifiably identify the very low-G behaviour as typical of a MV 
state [Fig. ©]. 

The characteristic feature of the Griffiths phases and the QCR is the local competi- 
tion between the MV and SEG states. This manifests itself as the co-existing clusters of 
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/-orbitals with mixed- and integer- valence. The integer-valence clusters are of greatest inter- 
est: these correspond to the dressing of a c-electron by a "cloud" of near-empty /-orbitals. 
This coupling of the c- and /-electron densities originates from the forward-scattering, and 
constitutes an electronic polaron. Electronic polarons form only at intermediate-coupling: 
at weak-coupling the screening clouds are heavily suppressed by the resonant (MV) scat- 
tering with the /-orbitals whereas at strong coupling the clouds merge to form the SEG 
phase. Since the bosonization method treats the forward-scattering exactly, our discovery 
of electronic polarons within a crossover regime between mixed- and integer-valence phases 
rigorously confirms the scenario proposed by Liu and Ho.— 

Note that in the Griffiths phases which flank the sides of the QCR the length scale £ 
differs from the conduction electron mean free path which gives rise to competing time 
scales: slow motion of the electronic polarons and fast motion of the conduction electrons. 
The different dynamics of the two types of particles provides a close analogy to a two fluid 
scenario. Since the polarons are randomly distributed across the lattice, these states can 
be viewed as intrinsic inhomogeneities involving charge fluctuations and short-range charge 
correlations. This resembles very closely the spin polaron liquid found in the ID Kondo 
lattice model.— 

Summarizing our results, we have derived an effective model for the occupation of the 
localized orbitals in the ID QFKM. This model predicts a crossover from a MV state at 
weak c — f Coulomb repulsion G to the integer- valence SEG state with increasing coupling 
strength. At intermediate values of the coupling the system exhibits an electronic Griffiths 
phase. The effective Hamiltonian gives a detailed accounting of the /-electron physics across 
the phase diagram, confirming the importance of electronic polarons in the QFKM. An 
extensive discussion of this system is to follow in a later paper—. 
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FIG. 1: Ground-state phase diagram for the QFKM below half- filling. The quantum critical regime 
(QCR) is found on either side of the solid line (5 = 0. Along the line n§ = 0.5 the QCR is replaced 
by a quantum critical point (QCP). The dotted lines denote the approximate "boundary" of the 
QCR with the Griffiths phases Gl and G2. The dashed lines separate these Griffiths phases from 
the segregated (SEG) and mixed-valence (MV) phases. Electronic polarons are present in both 
Griffiths phases and the QCR. 
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